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ABSTRACT 

We present a study, based on archival XMM-Newton observations, of the extended X- 
ray emission associated with the inner disk of M33. After the exclusion of point sources 
with Lx> 2 x 10 35 erg s _1 (0.3-6 keV), we investigate both the morphology and spectrum 
of the residual X-ray emission, comprising the integrated signal from unresolved discrete 
sources and diffuse components. This residual emission has a soft X-ray spectrum which 
can be fitted with a two-temperature thermal model, with kT w 0.2 keV and k 0.6 keV, 
the cooler component accounting for the bulk of the luminosity. There is evidence that the 
X-ray emitting plasma has a subsolar metal abundance. The soft X-ray surface brightness 
distribution shows a strong correlation with FUV emission and since the latter serves as a 
tracer of the inner spiral arms of M33, this is indicative of a close connection between recent 
star-formation activity and the production of soft X-rays. Within 3.5 kpc of the nucleus of 
M33, the soft X-ray and FUV surface brightness distributions exhibit similar radial profiles. 
The implication is that the ratio of the soft X-ray luminosity (0.3-2.0 keV) per unit disk area 
to the star formation rate (SFR) per unit disk area remains fairly constant within this inner 
disk region. We derive a value for this ratio of 1 — 1.5 X 10 39 erg s _1 ( Moyr -1 ) -1 , which 
is towards the top of the range of similar estimates for several other nearby face-on spiral 
galaxies {e.g. M51, M83). In the same region, the ratio of soft X-ray luminosity to stellar mass 
(the latter derived from K-band photometry) is 4 x 10 28 erg s _1 M _1 , a factor of 5-10 higher 
than is typical of dwarf elliptical galaxies (e.g. M32, NGC3379), suggesting that 10-20% of 
the unresolved emission seen in M33 may originate in its old stellar population. The remainder 
of the soft X-ray emission is found to be equally split between two spatial components, one 
which closely traces the spiral arms of the galaxy and the other more smoothly distributed 
across the inner disk of M33. The former must represent a highly clumped low-filling factor 
component linked to sites of recent or ongoing star formation, encompassing Hll regions, X- 
ray source complexes and X-ray superbubbles, whereas the distribution of the latter gives few 
clues as to its exact origin. 
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1 INTRODUCTION 

On the basis of current observations, the X-ray emission of spiral 
galaxies may be separated into several components. A very sig- 
nificant contribution generally arises from a set of resolved point 
sources which correspond to the most luminous examples of the 
galactic X-ray binary population, encompassing both high-mass 
(HMXB) and low-mass (LMXB) systems. To this must be added 
the integrated emission from large numbers of lower-luminosity 
sources, including supernova remnants (SNRs), more quiescent 
forms of X-ray binary, cataclysmic variables, and coronally active 
stars. Finally we need to include the thermal X-ray signal ema- 
nating from concentrations of truly diffuse million-K gas resident 
both in the galactic disk and potentially extending into the galactic 



halo. The extent to which these components can be distinguished 
from each other depends on the orientation and distance of the tar- 
get galaxy and, of course, on the spatial resolution and threshold 
sensitivity afforded by the observation. 

Early X-ray studies based on Einstein observations focussed 
on the brightest point sources seen in a number of nearby galax- 
ies, whilst also finding un derlying extended emission in complex 
structures (Fabbiano 1989). The superior spatial resolution and soft 
response of Rosat further progressed the field, for example allow- 
ing the extent of the apparently diffuse emission to be traced both 
in and out of the plane of the disk for face-on and edge- on systems 
respectively {e.g. iRead et alj[l997l ; iDahlem et all \l 998). More re- 
cent investigations, which exploit the enhanced collecting area, spa- 
tial resolution and spectral sensitivity of XMM-Newton and Chan- 
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dra, have probed the point-source populations in galaxies to much 



fainter thresholds than pre viously attainable (e.g. IStrickland et all 
l2004al : IColbert et alj2004h . The X-ray luminosity functions (XLF) 
typical of the HMXB and LMXB populations seen in spiral galax- 
ies have thereby been determined jGrimm et al]20 05; Kilear det al.l 
l2005l ; lFabbianol2006h . By extrapolation of the XLF to fluxes below 
the threshold at which individual sources can be resolved, the total 
aggregat ed luminosity of a p articular class of source may also be 
inferred dSazonov et al.ll200q) 

In two previous pa pers dWarwick et ai]|2007| . hereafter W07; 
lOwen & Warwickl2009l hereafter OW09), we studied the morphol- 
ogy of the underlying unresolved soft X-ray emission in a sample 
of six nearby face-on spiral galaxies and made a comparison with 
other wavelength, particularly FUV, measurements. These studies 
confirmed the close connection between recent star-formation and 
the production of soft X-rays. We also investigated how the soft 
X-ray luminosity to star-formation rate (SFR) ratio varied radially 
within the individual galaxies as well as from galaxy to galaxy, and 
found some evidence for an enhanced efficiency in the production 
of soft X-rays in regions of high SFR density. The link between 
soft X-rays and star formation can be further explored by studying 
galaxies closer to us, in which the majority of HMXBs and LMXBs 
can be readily resolved, so as to provide a relatively clear perspec- 
tive of the underlying soft X-ray emission. It is in this context that 
we focus in the present paper on M33. 

M33 is an Sc spira l gala xy with an inclination of 
56° dZaritskv. Elston & ffiiil 1 19891) at a distance of 795 kpc 
(van den Burgh 1991). As the the third largest galaxy in the Lo- 
cal Group, the relative proximity of M33 permits detailed study 
of its discrete X-ray source population and enables the separa- 
tion of these sources from any residual emission in the disk to a 
relatively faint threshold. The relatively low G alactic foreground 
N H (7.5 x 10 20 cm" 2 ; iKalberla et aljFiool) in the direction of 
M33, the low to moderate inclination of the galaxy, and its rel- 
atively high SFR in comparison to other nearby systems (0.3- 
0.7 M vr~ 1 . Iffippelein et aT]|2003h make M33 an ideal candidate 
for the present study. 

Early Einstein observations of M33 revealed a total of 17 
point sources, including the bright ULX M33 X-8 JEong et al.l 
198ll;lMarkert & Rallislll983| ; iTrinchieri. Fabbiano & Pereslll988l) . 



Trinchieri, Fabbiano & Peresl dl988h also found evidence for the 



presence of a soft diffuse component in the plane of the galactic 
disk. ROSAT observations dSchulman & Bregmanll995l : lLong et al.l 
1 1996b expanded the known populati on of X-ray point source s in the 
direction of M33 to a total of 184 dHaberl & Pietschl l200li). More 
recently, using XMM-Newton data, a total of 350-400 X-ray point 
sources to a limiting X-ra y luminosity of 10 35 erg; s -1 h ave been 
identified and categ o rized ( Foschini et al. 2004; Pietsch et 
iGrimm et al] 120051 ; iMisanovic et al.l 120061 ; iGrimm et ail 
These results h ave been complemented by recent studies o f M33 
using Chandra dPlucinskv et alj|2008l;lwilliams et alj|2008l) which 
have increased the source statistics particularly in the confused in- 
ner regions of the galaxy. Most of the emphasis to date has been 
on the X-ray emission from po int sources, bu t lHaberl & Pietschl 
d200ll) . IPietsch et alj l2004h and lTullman et alj d2008h have noted 
the presence of seemingly diffuse structures along the spiral arms. 
It is this latter component which we examine in this paper. 

Here, we focus on the spectral and spatial properties of the 
unresolved X-ray emission from the inner disk of M33, deduced 
from archival XMM-Newton observations. In fj2] we describe the 
set of observations used to construct soft-band X-ray images and 
outline the methods used in the subsequent data analysis. In Sj3] we 



briefly examine the properties of the resolved discrete point source 
population and the spatial distribution of the residual soft X-ray 
emission, which remains after the contribution of these sources is 
excluded. We go on to compare the soft X-ray morphology with 
GALEX FUV measurements and 2MASS K band data. We follow 
this with spectral analysis of the bright point source population and 
unresolved residual emission ([Q. We examine the spatial extent 
of the X-ray emission in comparison with star formation data and 
the underlying mass distribution (Sj5} for M33, and compare these 
relationships to those found in other galaxies. Finally we discuss 
the implication of our results (jj6]l and summarize our conclusions 
(0- 



2 OBSERVATIONS AND DATA REDUCTION 

In two previous papers (W07, OW09) we reported the results of an 
XMM-Newton study of six nearby face-on galaxies, namely M101, 
M83, M51, M74, NGC300 and NGC3184. Here we use similar 
methods to incorporate M33 into our sample. Since the 25 diam- 
eter of M33 is 70', which is more than a factor two larger than 
the 30' field of view of the European Photon Imaging Camera 
(EPIC) instrument on XMM-Newton, multiple-pointings are needed 
to give g ood coverage of the inner disk of M33. The survey con- 
ducted by Pietsch et al. (2004) matches this requirement and in the 
present work we have accessed the revelant datasets via the XMM- 
Newton public archive. Details of the observations which comprise 
the M33 survey are summarised in TableQ] 

Data reduction was based on SAS v8.0. The datasets were 
screened for periods of high background through the creation of 
full-field 10-15 keV lightcurves. MOS data were excluded when 
the 10-15 keV count-rate in a 100 s bin exceeded 0.2 ct s _1 , 
whereas pn data were excluded when the count rate in the same 
waveband exceeded 2 ct s . The resulting exposure times after 
filtering the observations, as reported in TableQ] range from 1.5 to 
13 ks in the pn camera, with typically somewhat greater exposure 
in the individual MOS cameras. Several of the pn datasets which 
include the bright ULX near the nucleus of M33 are affected by 
Out Of Time (OoT) events. The SAS task epchain was used to 
create "OoT event lists", which simulate the spatial distribution of 
OoT events in the observations. From these event lists, images of 
the distribution of OoT events were produced. These images were 
subtracted from the raw images produced from the original event 
lists to produce images cleaned of OoT events, which could then 
be used for imaging analysis. Additionally, the target observations 
were examined to look for evidence of solar win d charge exchange 
(SWC X) emission, using the methods detailed in lCarter & Sembavl 
(2008). One observation was found to be highly contaminated (ob- 
sid:0141980201) and was thus excluded from the dataset. 

In OW09 it was shown how an appropriate spatial mask could 
be used to excise the bright source population from the residual 
emission, thereby greatly reducing the "contamination" due to the 
former. Here we employ a slightly modified version of this image 
processing technique, necessitated by the need to mosaic the data 
sets from the M33 survey observations. 



2.1 Image Construction 

The methods we have employed to produce the X-ray images of 
M33 mirror the procedures described in OW09. For consistency, 
we use the same soft, medium and hard energy bands (0.3-1 keV, 
1-2 keV and 2-6 keV respectively), although in the present paper 



Soft X-ray emission from the inner disk ofM33 3 




Figure 1. Top-left panel: Adaptively smoothed version of the XMM-Newton pn+MOS image of M33 in the soft (0.3- 1 keV) band. Top-right panel: The same 
image in the medium (1-2 keV) band. Middle-left panel: Simulated image of the bright sources in M33, with the source mask contours overlaid. The bright 
central region is dominated by the ULX M33 X-8. Middle-right panel: The residual emission in M33 in the 0.3-1 keV band obtained by subtracting the 
bright source model and applying the spatial mask. Bottom-left panel: The GALEX FUV (X e f r ~ 1528A) image on the same spatial scale as the X-ray data. 
Bottom-right panel: 2MASS K band image on the same spatial scale as the X-ray data. The ellipse shown in each image represents the area of the galaxy used 
for X-ray analysis (see text for details). All the images are displayed with logarithmic amplitude scaling and are 40'on a side. 
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Table 1. Details of the XMM-Newton observations of M33. 



Galaxy 


Observation ID 


Start Date 


Filter a 


Target co-ordinates 


Useful exposure (ks) 






(yyyy-mm-dd) 


pn/MOSl/MOS2 


RA (J2000) 


Dec (J2000) 


pn 


MOS 1+2 


M33 


0102640101 c 


2000-08-04 


Ul-I- 


01 33 m 51.0 s 


+30° 39'37" 


7.1 


" 




0102640201 


2000-08-04 


M/M/M 


01 h 34 m 40.0 s 


+30° 57'48" 


11.8 


31.5 




0102640301 


2000-08-07 


M/M/Tn 


01 33 32.0 


+30 52 13 


3.6 


9.5 




0102640401 


2000-08-02 


Tk/Tk/Tk 


01 h 32 m 51.0 s 


+30° 36'49" 


9.1 


23.1 




0102640501 


2001-07-05 


M/M/M 


01 h 33 m 02.0 s 


+30° 21'24" 


9.2 


23.1 




0102640601 


2001-07-05 


M/M/M 


0l' l 34 m 08.0 s 


+30° 46'06" 


4.5 


11.9 




0102640701 


2001-07-05 


M/M/M 


01 34 10.0 


+30 27 00 


6.9 


21.9 




0102640801 


2001-07-07 


-/M/M 


0l' l 34" l 51.0 s 


+30° 42'22" 




3.2 




0102640901 


2001-07-08 


M/M/M 


01 ft 34 m 04.0 s 


+30° 57'25" 


3.9 


11.2 




0102641001 


2001-07-08 


M/M/M 


0l' l 33 m 07.0 s 


+30° 45'02" 


1.5 


16.3 




0102641101 


2001-07-08 


M/M/M 


01 h 32 m 46 qs 


+30° 28' 19" 


8.0 


21.0 




0102641201 


2000-08-02 


Tk/Tk/Tk 


0l' l 33 m 38.0 s 


+30° 21'49" 


12.0 


7.2 




0102642001 


2001-08-15 


M/M/M 


0l' l 34 m 51.0 s 


+30° 42'22" 


8.8 


22.3 




0102642101 


2002-01-25 


M/M/M 


01 h 34 m 34.0 s 


+30° 34'11" 


10.0 


24.3 




0102642201 


2002-01-25 


M/M/M 


01 h 34 m 56 qs 


+30° 50'52" 


11.5 


27.3 




0102642301 c 


2002-01-27 


M/M/M 


0l' l 33 m 33.0 s 


+30° 33'07" 


9.9 


24.1 




0141980101 


2003-07-11 


M/M/M 


0l' l 33 m 07.0 s 


+30° 45'02" 


6.2 


13.2 




0141980301 


2003-07-25 


-/M/M 


01 h 34 m 08 qs 


+30° 46'06" 




1.2 




0141980501 


2003-01-22 


M/M/M 


0l' l 33 m 51.0 s 


+30° 39'37" 


1.9 


17.0 




0141980601 


2003-01-23 


M/M/M 


0l' l 32 m 51.0 s 


+30° 36'49" 


11.0 


25.9 




0141980701 


2003-01-24 


M/M/M 


0l' l 33 m 38.0 s 


+30° 21'49" 


4.4 


11.4 




0141980801 c 


2003-02-12 


M/M/M 


0l' l 33 m 51.0 s 


+30° 39'37" 


7.8 


19.8 



a - Tn = thin filter, M = medium filter, Tk = thick filter 

- Assumed position of the galactic nucleus. 
c - Observations used for spectral analysis. 



emphasis is placed on the soft band, which contains the majority of 
the diffuse signal. Images and exposure maps were extracted for all 
three bands and cameras for each observation. A constant particle 
rate, estimated from the corners of the detector not exposed to the 
sky, was subtracted from each image and the set of images for a 
given bandpass were co-added using the method detailed in OW09, 
with a position offset dependent on the pointing direction of the 
observation. As the majority of observations were conducted using 
the medium filter, factors were derived to convert observations with 
the thin and thick filters to the same relative intensity scale. The ex- 
posure maps were likewise mosaiced to produce a single exposure 
map for the extended field. Finally, flat-fielding was achieved by di- 
viding the mosaiced image for each band by the relevant exposure 
map. 

At this stage it proved necessary to apply an iterative adjust- 
ment to the particle background estimates, so as to create a visu- 
ally flat image. This was achieved by imposing a mask on the flat- 
fielded images consisting of a point-source mask (see S]2.2| l cou- 
pled with an elliptical mask extending to a major axis radius of 
15' from the nucleus of the galaxy. The ellipticity of the latter was 
fixed at the major/minor ax is ratio reported in the RC3 catalogue 
dde Vaucouleurs et al]|l99lh for M33 , that is 1.7:1, with the posi - 
tion angle of the major axis set to 22° dCorbelri & Walterboj2007h . 
We next imposed the requirement that the level of the X-ray sky 
background should be constant in the flat-fielded images outside 
of the masked region (i.e. we assume that there are no significant 
gradients in soft X-ray background across the full field of the sur- 
vey and also that there is negligible contribution from M33 in this 
outer region). This did in fact appear to be justified based on visual 
inspection of the raw images. Next, the average level of the sky 
background in the un-masked region of the image was calculated. 
Using this average value, the particle background rates for each 



observation then were revised, so as to force the sky background 
outside of the masked region in each observation to this average. 
The resulting mosaiced and flat-fielded images of M33 in the soft 
and medium bands are shown in Fig. [I] 

We note that given the above, radial analysis of the surface 
brightness distributions contained within the M33 images is only 
justified for the region extending from the nucleus to a major-axis 
radius of 15'. This equates to 42% of the 25 region and to a radial 
linear extent of 3.5 kpc at the distance of M33. We describe this 
region hereafter as the inner disk of M33. 

2.2 Spatial masking of bright sources 

The catalogue of sources used to produce a bright- source spatia l 
mask for M33 was taken from the work of lMisanovic et al, I feOOfjh . 
who reported 350 X-ray sources across the 25 disk of M33 above 
a luminosity of 2 x 10 erg s . Ninety-two of these sources lie 
within the elliptical region defined in £]2.1l (see Table [2] for sum- 
mary details). Following the methods described in OW09, a model 
"bright-source image" was created for M33. A surface brightness 
cut was applied to this image at a level of 0.07 pn+MOSl+MOS2 
ct ks - pixel -1 to produce the source mask (see Fig[T). We esti- 
mate the "spillover" fraction of this mask, i.e. the fraction of the 
bright-source signal not contained within the masked region, to be 
« 4%. This mask was used to divide the inner disk region into 
two components, namely a "bright-source region" and a "residual- 
emission region" (using the terminology of OW09), which could 
then be subject to appropriate spatial and spectral analysis^. In the 

1 The source mask is square in the region of the very bright central source 
due to clipping of the specified surface brightness contour at the edge of the 
4' X 4' sub-image used to represent the point spread function. The leakage 
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image analysis, low-level contamination of the residual-emission 
signal by the spillover from the bright sources was further sup- 
pressed by subtracting the simulated image from the correspond- 
ing pn+MOS image and reimposing the source mask. In contrast, 
for spectral analysis, the spillover was accounted for by including 
a bright-source contribution in the spectral fit of the residual emis- 
sion. 

2.3 Spectral Extraction 

The soft-band image in Fig. [T] demonstrates the existence of an ex- 
tended emission component in addition to the population of bright 
point sources. On the basis of the approach described earlier, we ex- 
tracted the integrated pn spectrum of both the bright-source region 
(bounded by the spatial mask) and the residual-emission region 
(corresponding to the full X-ray extraction region less the source- 
masked area). This process was carried out for the three observa- 
tions indicated in Table [T] which were chosen as they encompass 
the central region of the galaxy with reasonable observation times. 
A narrow strip of one CCD contaminated by Out of Time (OoT) 
events from the ULX was excised in each observation. The SAS 
tools arfgen and rmfgen were used to produce appropriate Auxil- 
iary Response File (ARF) and Response Matrix File (RMF) files for 
the source and residual galaxy regions, and the counts recorded in 
adjacent (raw) spectral channels were summed to give a minimum 
of 20 counts per spectral bin in the final set of spectra. 

The large extent of the bright-source and residual-emission re- 
gions makes the process of determining appropriate background 
spectra more complex than for most XMM-Newton applications. In 
OW09, we used appropriately scaled spectra from an annulus sur- 
rounding the defined galaxy region and the corner regions of the 
detector to approximate the background. This process is not viable 
for M33 as the central galaxy field covers too large a fraction of the 
EPIC pn field of view, thus making the extraction of an adequate 
background region impossible. We therefore used a combination of 
"blank-sky" fields extracted from a region of sky close to M33 (to 
minimize the difference in the sky X-ray background) and "filter- 
wheel closed" data to produce a background spectrum. Using the 
SAS tool skycast, a co-added blank-sky pn image was rotated to the 
same attitude as each observation, and the spatial mask imposed for 
the bright source and residual regions. Spectra were then extracted 
from these areas and scaled to the same exposure time as the ac- 
tual observation. The same process was followed for filter-wheel 
closed data to produce a particle background spectrum. The source 
and blank-sky spectra were then compared in the 8-12 keV band, 
where the signal is dominated by the particle background. The dif- 
ference in this signal was compensated for by addition of a relevant 
fraction of the filter-wheel closed spectrum, the effect of which was 
to ensure the particle background level for source and background 
spectra was the same. 



3 PROPERTIES OF THE GALACTIC X-RAY EMISSION 

3.1 The contribution of luminous point sources 

The relative proximity of M33 means that point sources in this 
galaxy can be resolved (and excluded) down to a luminosity thresh- 
old of 2 x 10 35 erg s _1 (0.3-6 keV), which is ~ 100 times deeper 

of flux from the bright central source is calculated with respect to the mask 
actually employed. 



than was the case for the majority galaxies in our earlier study 
(OW09). A significant fraction of X-ray luminosity deriving from 
HMXBs and LMXBs in M33 can thereby be excluded. Using a 
conversion from count-rate to luminosity derived from the spectral 
fit of the bright point source population in the galaxy, we measure 
the summed X-ray luminosity (0.3-6 keV) of the bright sources 
contained within the source mask to be 1.6 x 10 39 ergs -1 , to 
which the bright ULX, M33 X-8 makes a dominant contribution. 
We determine the X-ray luminosity of this individual source (av- 
eraged across the set of M33 survey observations) to be 1.2 x 
10 39 erg s _1 (0. 3-6 keV), bro adly consist ent with the measure - 
ments reported bv l Grimm et alj d2005h and lMisanovic et all (2006), 
who report a luminosity of 8.3 x 10 38 erg s _1 (0.2^1.5 keV). This 
source is detected within 0.6" of the optical nucleus of the galaxy 
dDubus & Rutledgd 12002). but short-term variability in its X-ray 
e mission suggests that it is not the galactic nucleus itself. Analysis 
bv lFoschini et al.1 ( 120041) suggests that the object is a black hole of 
mass «10 Mq accreting at a super-Eddington rate. 

Table [3] summarizes the distribution across the soft, medium 
and hard energy bands of the integrated luminosity contained 
within the bright-source region (with the conversion from count 
rate to flux and hence luminosity estimated from the best-fitting 
spectral models reported in §4.2\ . The X-ray luminosity pertain- 
ing to the residual-emission component is similarly reported. The 
quoted Lx figures are corrected for foreground galactic absorption. 
Correction factors were also applied for spillover of source counts 
into the residual-emission region and for the underlying extended- 
emission component contained within the source mask. 

It is possible to estimate the integrated X-ray luminos- 
ity of point sources below our luminosity threshold, assuming 
such sources have similar characteristics to the resolved source 
population. For thi s, we use the Chandra M33 observations of 
iGrimm et alj 20051). who s t udy a similar area of the galaxy as con- 
sidered here. Grimm et al.1 d2005l) derive an XLF with a slope of 
-(0.74—0.78), which we use to estimate the total integrated Lx for 
sources between 2 x 10 35 erg s _1 and 1 x 10 34 erg s _1 . The re- 
sults are given in Table [3] and show that unresolved X-ray binaries 
do not contribute significantly to the residual X-ray emission ob- 
served in the soft band, but may form a significant contribution to 
the emission observed above 1 keV. This implies that the residual- 
emission component largely arises as the integrated emission of 
lower-luminosity source populations such as cataclysmic variables 
and active binaries plus truly diffuse emission associated with the 
inner disk of M33, presumably energized by the collective effect of 
supernovae explosions and stellar winds. 

3.2 Morphology of the soft X-ray emission 

M33 contains a two-armed grand-design spiral structure with spi- 
ral arm s extending to a radiu s of 10' as observed in the near in- 
frared (Regan & Vogel 1994), along with a bar structure in the 
inner 1.5'. Far infrared and radio observations confirm this struc- 
ture, with the southern spiral arm the most promin ent emission re- 
gion in the galaxy across al l of these wavebands ( Hippele in et al.1 
2003; Tabatabaei et al. 2007). The soft X -ray image of M33 (Fig.[T]( 
clearly shows substantial, apparently diffuse emission distributed 
across the inner disk of the galaxy, with the overlaying pattern 
of the spiral structure particularly pronounced in the region of the 
southern arm. Narrow X-ray structures are similarly observed along 
the northern spiral arm, although they are noticeably less luminous. 
The soft X-ray emission can be traced out to a radius of 15' beyond 
which the signal drops below the level of the X-ray background. 
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Table 2. The parameters of the region investigated in this paper. 



Galaxy 


X-ray extent" 
(') 


Major-axis PA 

(°) 


Threshold L x 6 

(10 35 ergs" 1 ) 


Number in 
Source List 


Number of high Lx sources 
(L x > 5 X 10 38 ergs" 1 ) 


M33 


30.0/17.6 


22 


2.0 


92 


1 



a - The major/minor axis diameters of the "X-ray extraction region". 
b - Nominal Lx threshold applied in the 0.3-6 keV band in defining the bright-source sample. 



Table 3. Breakdown of the total X-ray luminosity of M33. 



Galaxy Spillover/Area Component Lx(10 38 erg s 1 ) 

Factors (%) (0.3-1 keV) (1-2 keV) (2-6 keV) (0.3-6 keV) 



M33 


4/6 


Bright sources" 


5.6 


4.0 


6.2 


15.8 






Unresolved sources 6 


[0.06] 


[0.04] 


[0.07] 








Residual emission 


1.1 


0.1 




1.2 






Total measured 


6.7 


4.1 


6.2 


17.0 



a - Down to a threshold Lx= 2 X 10 ib erg s . 
b - Extrapolated to the lower limit Lx= 1 x 10 34 erg s" 



Fig. [T] shows a comparison of the soft (0.3-1 keV) X-ray 
morphology of M33 with corresponding FUV (\ e ff ~ 1528A) 
and near-infrared K-band measurements. The FUV data are from 



GALEX jGil de Paz et alj |2007[) and the K-band data from the 
2MASS survey JSkrutskieet ail 2006; mosaic images provided by 
T. Jarrett (IPAC)). The GALEX and K-band images were resampled 
to the same spatial scale as the XMM-Newton images, and were 
then lightly smoothed with a gaussian mask with a ~ 4" (roughly 
matching the XMM point-spread function). The FUV emission 
traces the locations of recent star formation, which in turn map 
out an underlying spiral structure. In contrast the K-band emission 
shows a relatively smooth azimuthal distribution coupled with a 
sharp decrease in the surface brightness with increasing galacto- 
centric radius. The soft X-ray morphology appears to share some 
of the characteristics of both FUV and K-band images, and the as- 
sociation of the X-ray emission with each of these components is 
explored more quantitatively in Section|5] 



4 SPECTRAL ANALYSIS 

The methodology outlined in i]2.3l was employed to extract bright- 
source and residual-emission spectra for the three observations 
identified in Table [T] Only EPIC pn data were used in the spec- 
tral fitting, since the pn camera has superior sensitivity to the MOS 
cameras in the soft band. The bright-source spectrum was fitted 
separately for each observation to allow for possible spectral vari- 
ations between observations (and also to take account of the fact 
that some of the sources encompassed by one observation might 
fall outside the field of view in another), whereas the three residual- 
galaxy spectra were fitted simultaneously. The spectral fitting was 
carried out using the software package XSPEC version 12.5. 



4.1 Spectra of the bright-source region 

As noted earlier, the X-ray luminosity of the bright-source pop- 
ulation of M33 is dominated by M33 X-8, the ULX close to the 
nucleus of the galaxy. It follows that our bright-source spectra will 



also be dominated by this source. We conducted the spectral fit- 
ting by first extracting and then modelling the spectrum of the 
ULX alone using the combination of a multicolour black-body disk 
component and a pow er-law continuum, bo th subject to some in- 
trinsic absorption (see iFoschini et all 120041) . The spectrum of the 
ULX appears to change marginally between observations, reflect- 
ing the long-term spectral variability of this source. Next we mod- 
elled the bright-source spectra extracted from the full masked re- 
gion by combining an unabsorbed power-law component with the 
ULX contribution, the former representing the emission of the other 
bright sources (mainly LMXBs and HMXBs) in M33. Fig. |2] (top 
panel) illustrates the bright-source spectrum and derived composite 
best-fit spectral form for one of the observations. The separate con- 
tributions of M33 X-8 and the other bright sources to the composite 
bright-source spectrum are listed in Table|4]for each observation. 

We estimate that ~ 3% of the flux from the ULX leaks be- 
yond the source mask into the residual galaxy region, whereas the 
corresponding value for the other bright sources is ~ 9% (the dif- 
ference reflecting the use of a fixed surface brightness cut to define 
the source mask). These percentages were used when modelling the 
bright-source spillover into the residual-emission region. 



4.2 Spectra of the residual-emission region 

The three spectral datasets which were selected as representative 
of the residual emission are shown in Fig. [2] As a preliminary step 
these datasets were checked for short-term variability at ~ 0.6 keV 
indicative of the presence of OVII geocoronal SWCX emission. No 
evidence for such contamination was found and in fact the three ob- 
servations give a consistent measure of the strength of the OVII line 
when interpreted as thermal emission from M33 (see below). Nev- 
ertheless we cannot completely rule out the presence of some resid- 
ual OVII emission associated with heliospheric SWCX, for which 
the variability timescale would be relatively long. 

The residual-emission spectra were modelled with a com- 
bination of two thermal plasma (Mekal) components both sub- 
ject to line-of-sigh t absorption in the Galactic foreground, Nh 1 ^ 
7.5 x 10 20 cm -2 tealberla et al]|2005h . Initially the two thermal 
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Figure 2. The EPIC spectra for M33. Top panel: The spectrum of the bright-source regions from one observation. Bottom panel: The simultaneously fitted 
spectra of the residual-galaxy region. In all cases the solid line corresponds to the best-fit spectral model (see text). The \ 2 residuals with respect to the 
best-fitting model are also shown in both panels. 



components were constrained to solar metal abundances. A "cool" 
thermal component at kT « 0.2 keV, together with a "hot" ther- 
mal component at kT ~ 0.6 keV provided a reasonable fit to the 
spectral data (x^=l-05 for the fit to all three spectra). Allowing the 
metallicity of the two thermal components to vary independently 
resulted in a constrained fit only for the cool component. With the 
abundances of the two components tied, the best fit was obtained 
for a metallicity relative to solar of 10±3%. The resulting improve- 
ment in the \ 2 (X?=0-99 for the spectra set as a whole) was highly 
significant as measured by the F-test. Details of these spectral fit- 
ting results are summarised in Tableland the resulting best-fitting 
model spectrum (with subs olar abundances) is illustrat ed in Fig. [2] 
Several authors (e.g. IStrickland & Stevensl|l998i ; W07) have 
noted that when attempting to fit low-resolution X-ray spectra per- 
taining to complex multi-temperature plasmas with over-simplified 
one- or two-temperature models, then a requirement for strongly 
sub-solar abundances is often the outcome, alas an erroneous one. 
Notwithstanding this cautionary note, in the present case, the in- 
troduction of a sub-solar abundance does improve the spectral fit 



and more specifically removes a soft excess otherwise apparent in 
the residual-galaxy spectrum below 0.5 keV. This is evidence, al- 
beit tentative, that the low-metal abundance may not necessarily be 
a fitting artifact in this case. Interestingly the adoption of sub-solar 
abundances would also compensate for the apparent soft X-ray ex- 
cesses found previously in M51 and M83 (OW09), although for 
these sources the signal to noise ratio was insufficient to give use- 
ful individual constraints. In the case of M33, the subsolar metal- 
licity inferred from the soft X-ray spectroscopy is also in line with 
studies of its HII regions, for whi ch O/H abundances a factor of 2-3 
below solar have been reported ICrockett et al .I I 2006). 

The X-ray luminosity of the residual-emission component in 
M33 is 1.2 x 10 38 erg s" 1 (0.3-2 keV). This is approximately 10 
times lower than the value inferred for the majority of the galaxies 
studied by OW09 - consistent with M33 being a relatively low- 
mass system in a satellite orbit around M3 1 . The measurement of 
a thermal spectrum which may be characterised by temperatures of 
Rs 0.2 and « 0.6 keV is consiste nt with many previous results for 
normal and starburst galaxies (e.g. lEhle et al.ll998l : lFraternali et all 
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Table 4. Parameters of the best-fitting models to the spectra of the bright-source and residual-galaxy regions. 



Observation 



Region 



Intrinsic JVh 

™-2 



Power-law 
Index 
Normalization 



Disc BB 
keV 
Normalization 



Cool MEKAL 

keV 
Normalization 



Hot MEKAL 

keV 
Normalization 



Goodness CookHot 
of Fit Flux Ratio 

X 2 /dof (0.3-2 keV) 



0102640101 



Bright Source: 
M33 X-8 



9.3 x I0 20 2.61±0.04 1.16±0.01 
1.99 x 10~ 3 0.318 



Bright Source: 
Other 



1.70±0.06 

5.48 X 10- 4 



0102642301 



Composite Fit 

Bright Source: 
M33 X-8 



3.5 x 10 21 2.90±0.08 
3.06 x 10~ 3 



676/703 



1.23±0.03 
0.150 



Bright Source: 
Other 



2.04±0.02 
1.90 x 10" 3 



0141980801 



Composite Fit 

Bright Source: 
M33 X-8 



8.9 x 10 20 



2.30±0.03 
2.83 x 10" 3 



720/723 



1.16±0.02 
0.143 



Bright Source: 
Other 



6.5 X 10 20 



1.86±0.07 

5.50 x 10~ 4 



Composite Fit 



709/672 



Combined 



Residual Galaxy 
Solar-abundance 



0.17±0.01 

8.80 x 10" 4 



0.59±0.03 
1.10 x 10" 4 



1337/1266 



4.5 



Combined 



Residual Galaxy 
Free-abundance" 



0.18±0.01 

6.35 X 10~ 3 



0.63±0.07 
6.01 x 10" 4 



1254/1265 



4.2 



■ Residual emission fit at 10±3% solar abundance. 



Table 5. Physical properties of the diffuse gas present in each galaxy. 



Galaxy Radius" Component Electron Density Thermal Energy Cooling Timescale 
kpc (keV) lO" 3 ^- 1 / 2 cm" 3 10 5i r] 1 / 2 erg lO 8 ?? 1 / 2 yr 



M33 



3.5 



0.2 
0.6 



4.1 
1.4 



4.4 
4.8 



17 

61 



■ Assumed radius of a putative shallow halo component (see text) 



l2002l ; ISoria & Wul2003l : lKuntz et alj20oS OW09). In M33 the rel- 
ative contribution of the two thermal components in the 0.3-2 keV 
band is 4.2:1, with the cool component dominant. M33 is therefore 
spectrally similar to M74, M101 and NGC3184, as opposed to M51 
and M83, where there is a more equal weighting between the two 
thermal components. 

We can derive the physical properties of the diffuse gas from 
the derived spectral parameters. If we assume that the majority of 
the residual emission we observe is truly diffuse and is contained 
within a cylindrical disk of major-axis radius 15'(;.e. a linear di- 
mension of 3.5 kpc) and half-width 0.5 kpc, we can estimate the 
electron density, the thermal energy and cooling timescale for each 
thermal component. The results are summarised in Table [5] The 
thermal energy contained in the two components is comparable, 



implying that these components they may be in rough pressure bal- 



5 SPATIAL ANALYSIS 

In an earlier study (W07), we found that when we compared the 
observed soft X-ray surface brightness distribution of M101 with 
corresponding datasets measured in several optical to ultraviolet 
wavebands, the best correlation was obtained with the U-band im- 
age. This was interpreted in terms of an X-ray signal comprised of 
two distinct spatial components, namely a clumpy thin-disk com- 
ponent which traces the spiral arms of the galaxy, and an underly- 
ing spatially-smooth component which contributes significantly to 
the central concentration of the soft X-ray emission. The soft X-ray 
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morphology of M33 (see §3) appears to bear some similarity to that 
of M101 in that it is possible to discern both the spiral arm struc- 
ture (which clearly dominates the corresponding FUV image) and 
a smoother underlying component (a template for which might be 
provided by the K-band image). 

To explore the above idea further, we have investigated the de- 
gree to which a "synthetic X-ray image" produced as a linear com- 
bination of the FUV and K-band images can be made to mimic the 
observed soft X-ray emission. To this end, the soft X-ray, FUV and 
K-band images shown in FigQ]were all spatially masked so as to si- 
multaneously suppress the bright X-ray point sources, several FUV 
regions with very high surface brightness and some bright fore- 
ground stars visible in the K-band image. The three images were 
then compressed to a 1' x 1' pixel scale. A synthetic image formed 
by summing scaled-versions of the FUV and K-band images was 
then compared to the corresponding soft X-ray image through the 
computation of a x 2 statistic (calculated on a pixel-by-pixel basis 
and summed over the set of pixels within the elliptical inner disk 
region defined earlier). By varying the relative contributions of the 
FUV and K-band images to the synthetic image we were thereby 
able to identify the best-fitting combination. The best fit was in fact 
obtained when the FUV template image contributed 40% of the to- 
tal signal (in the synthetic image) and the K-band equivalently 60% 
(see Fig[3](. This analysis confirms a direct association between re- 
cent star formation and soft X-ray emission in the inner galactic 
disk of M33. However, it also points to a substantial contribution 
from an underlying component correlated with the K-band light, 
which in turn traces the old stellar population of the galaxy and, to 
first order, the galactic stellar-mass distribution (see §6). 

In OW09, we derived the azimuthally-averaged radial profiles 
of the soft X-ray emission observed in six late-type spiral galaxies 
and used these measurements, in conjunction with published SFR 
data, to derive quantitative estimates of the ratio of the soft X-ray 
emission per unit disk area to the local star formation rate (SFR) 
per unit disk area. Here we determine the same information for 
M33 using the same methodology as OW09, except that since M33 
is of intermediate inclination, with a major/minor axis ratio of 1.7, 
in this case it is necessary to use elliptical annuli centred on the 
galactic nucleus in order to determine the radial profiles. Due to the 
presence of the ULX near the galactic nucleus of M33, we extracted 
data from the spatially-masked soft X-ray image between major- 
axis radii of 2' and 14', beyond which the soft X-ray emission falls 
to near the background level. The corresponding range of linear 
extent at the distance of M33 is 0.46 to 3.2 kpc. In the analysis 
w hich follows, we use the FUV and K-band radial profiles tabulated 
in lMunoz-Mateos et al] < l2007l) . 

The K-band and soft X-ray radial profiles (FigO top panel) 
show very similar rates of fall-off between 0.5 and 3 kpc; 
iMunoz-Mateos et al. U20071) quote a scale length of 1 .48 kpc for the 
K-band photometry. In contrast, the FUV radial profile although of 
comparable slope to the X-ray and K-band curves outside of 2 kpc, 
appears to flatten at smaller radii. Although, given the background 
uncertainties, we are unable to place quantitative limits on the level 
of the soft X-ray emission much beyond 3 kpc, visual inspection of 
the X-ray images suggest that the soft X-ray signal declines very 
sharply outside of the inner-disk region. This rapid switch-off is 
reminiscent of a similar effect evident in NGC300, M74, M5 1 and 
M101 (OW09). 

Fig. ID (second panel) shows a comparison of the soft X-ray 
luminosity in the 0.3-1 keV band (per unit di sk area) with the SFR 
(per u nit disk area) (the latter as tabulated in Munoz-Mateos et al. 
2007). More specifically the ratio of these two quantities has a value 



of ?s 1 x 10 39 erg s _1 ( M Q yr" 1 ) -1 , at a radius 0.5 kpc, with ev- 
idence for a modest decline (by a factor dex~0.2) across the sam- 
pled region. In order to make a direct comparison with the results 
for the galaxies studied in OW09, we need to apply two corrections 
as follows: (i) we convert the X-ray luminosities to a broader 0.3- 
2 keV band using the spectral models discussed in W.2I - this re- 
sults in an upward scaling of the X-ray/SFR ratio by 5%); (ii) we 
correct the X-ray luminosities to a common point source exclusion 
threshold of L x = 1 x 10 37 erg s _1 (0.3-6 keV) (in the case of M33 
this involves integrating the luminosity encompassed by its point- 
source XLF over the range Lx= 2 — 100 x 10 35 erg s _1 ) - this 
results in a further upward scaling of 50%. The resulting estimate 
of the X-ray/SFR ratio of 1 - 1.5 x 10 39 erg s" 1 ( M Q yr -1 ) -1 
lies towards the upper end of the range observed to date, at a very 
similar level to that pertaining to the inner disks of M5 1 and M83 
(cf. Fig.5, OW09, where for M33 the parameter log(SFR density) 
ranges from -8 to -7.50. 

The soft X-ray radial profile and the stellar-mass radial pro- 
file est imated from the K-ba nd photometry can similarly be com- 
pared. IMunoz-Mateos et aD §007) use K-band data to evalu- 
ate the stellar-mass radial profile of M33. We find that the soft 
X-ray/stellar-mass ratio remains relatively constant at w 4 x 
10 28 erg s" 1 M _1 , across the inn er disk of M33 (Fig. [3 third 
panel). Recent studies conducted by Revnivtsev et alj d2007l) and 
Revnivtsev, Vikhlinin & Sazono 3 J2007I) have explored the rela- 
tionship between unresolved X-ray emission and the stellar-mass 
distributions for the Galactic ridge and M32, in an attempt to de- 
rive the X-ray emissivity per unit stellar mass in systems with low 
SFR. Most likely this emissivity originates in point sources associ- 
ated with (intermediate) disk population objects such as active bina- 
ries and cataclysmic variables with Lx between 10 30 erg s _1 and 
10 34 erg s _1 , together with a contribution from fainter sources, 
the bulk of which will be coronally acti ve stars. From these an d 
other studies (e.g. NGC33 79 and M31 - iRevnivtsev et al.1 < |2008|) : 
iBogdan & Gilfanovl §008)), an X-ray to stellar mass ratio for qui- 
escent systems (which do not contain substantial hot gas) can be 
estimated to be 4 - 8 x 10 27 erg s _1 M© -1 (0.5-2 keV). This is 
a factor of 5-10 lower than we derive above for M33, implying 
that m 10 — 20% of the observed unresolved X-ray luminosity in 
M33 may be associated with the old stellar population. This level 
is consistent with that predicted from the extended XLF derived 
for LMXB, relatively quiescent LMXB, cataclysmic variables and 
other source types associated with the old population. 

As a final investigation we also considered the variation in X- 
ray spectral hardness as deduced from images produced in two sub- 
bands, namely the 0.3-0.8 keV and 0.8-1.2 keV bands. The plot of 
spectral hardness versus radius shown in (Fig. [4] bottom panel) in 
essence tracks the relative contribution of the two thermal spectral 
components identified earlier in the spectral fitting ( i]4.2t . There is 
a hint of spectral softening with increasing galactocentric radius, 
suggesting that the hotter thermal component makes a greater con- 
tribution to the total luminosity towards the centre of the galaxy. 



2 In the discussion section of OW09, the upper bound of the X-ray to SFR 
ratio was quoted as 5 X 10 39 erg s — 1 ( Mq yr — This included an 
additional scaling factor of 3.4 arising from the conve rsion from the SFR 
mass range employed by Munoz-Mateos et al. (2007) (0.1-100 Mq) to 
that used in|Mas-Hesse, Oti-Flo ranes & Cervind 120081) (2-120 Mq). 
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Figure 3. Left-hand panel: \ 2 contour plot for the match of the synthetic image derived from the FUV and K-band templates to the observed soft X-ray surface 
brightness distribution. The best fit is obtained when the split between the FUV and K-band components is approximately 40:60 (marked with the cross). The 
contours represent 68%, 90% and 99% \ 2 confidence intervals. Right-hand panel: The synthetic image which best matches the soft X-ray morphology. The 
displayed image has 4" pixels whereas the \ 2 tes t was performed using 1' binning (see text). The amplitude scaling is logarithmic. The ellipse has a major-axis 
radius of 15'and corresponds to the X-ray analysis region. 



6 DISCUSSION 

The strong correlation observed between soft X-ray and FUV 
emission in M33 confirms a direct linkage between star forma- 
tion and the production of soft X-rays in this galaxy. Accord- 
ing to current models (e.g. Cervino. Mas-Hesse & Kunthl |2002| : 



IMas-Hesse. Oti-Floranes & Cervinol l2008h . following a burst of 
star formation which gives rise to a massive star cluster, the FUV 
emission may be expected to peak on a timescale matching the life- 
time of the most massive stars (~ 1 — 3 x 10 6 years), whereas 
the X-ray signal reaches a maximum sometime later. In this sce- 
nario, diffuse soft X-rays are produced through the heating of bub- 
bles within the interstellar medium to temperatures of 10 6-7 K as 
a result of the mechanical energy input from the winds of mas- 
sive stars created in the starburst and the eventual destruction of 
such stars in supernovae. Individual SNR also contribute, primar- 
ily during the free-expansion and adiabatic phases of their evo- 
lution. Occasionally luminous high-mass X-ray binaries might be 
formed in such a cluster, although here we assume that such sources 
would be bright enough to be excluded as resolved point-sources. 
IMas-Hesse, Oti-Floranes & Cervinol [2008 find that following the 
onset of the starburst the diffuse soft X-ray luminosity increases 
rapidly over the first few Myr. For an instanteous burst, the X- 
ray luminosity may then plateau and eventually decline whereas, 
if the starburst activity is ongoing then the X-ray luminosity may 
gradually rise for up to ~ 30 Myr until the stellar formation 
an d destruction process reach an equi l ibrium . Using the models 
of IMas-Hesse. Oti-Floranes & Cervinol d2008h . we find that the X- 
ray/SFR ratio measured for the inner disk of M33 matches the pre- 
dictions at ~ 10 Myr (after the onset of an extended burst of star 
formation activity), assuming the efficiency of conversion of me- 
chanical energy from supernovae into X-rays is 1% (cf. O W09). 
Studies of the star format ion history of M33 ^Wilson et al.ll 19881 ; 
I Wilson & Matthews| [l995) have shown that several strongly emit- 
ting HII regions have undergone bursts of star formation on the 
timescale of 10 Myr, which indicates that this model is realistic. 



Through a quantitative comparison of the observed soft X-ray 
surface brightness distribution with a synthetic image, we found 
that roughly 40% of the X-ray signal had a spatial distribution sim- 
ilar to that seen in the FUV with the remaining 60% better match- 
ing the K-band distribution. In effect this represents a 40:60 split 
between a clumpy spiral-arm distribution and a smooth disk dis- 
tribution. However, in §6 we also discussed the fact that perhaps 
20% of the total observed soft X-ray emission might be directly 
associated with the old stellar population of the galaxy (for which 
the K-band light again serves as a tracer). If we subtract this frac- 
tion we are left with a roughly equal split between the two inferred 
spatial components of the soft X-ray emission. 

We are now in a position to add some refinements to our 
earlier analysis in which we modelled the X-ray emitting gas as 
filling a cylindrical disk of half-thickness 0.5 pc, extending to a 
galactocentric radius of 3.5 kpc. As the spiral arm regions oc- 
cupy at most ~ 20% of the inner galactic disk by area, we 
may model 40% of the emission (directly connected to the spi- 
ral arm regions) as filling 20% of the disk with a thickness of, 
say, 200 pc. With this adjustment, the cooling timescale of the 
gas is ~ 2 x lO 8 ^ 1 ' 2 yr, where r\ is the filling factor of the 
ga s. On th e basis of the dynamical analysis of M33 performed 
bv lPueraril Jl993h . narrow features initially distributed in a spi- 
ral arm pattern would be completely smeared out by differential 
galactic rotation on timescales of 2 — 3 x 10 8 yr. For such nar- 
row features to survive we need a much shorter cooling timescale 
(assuming radiative cooling is the dominant process) implying a 
filling factor much less than unity. Setting r\ ~ 10~ 3 gives a cool- 
ing timescale of 6 x 10 6 years, comparable t o the inferred lifetime 
of ho t bubbles observed in our own Galaxy (Egaer & AschenbachI 
1 19951 : Breitschwerdt, de Avillez & Baumgartner l2009h . As noted 
earlier, this is also roughly the timescale on which the soft X- 
ray production maximises following the onset of a starformation 
episode. Given the rotational period and pattern speed of the galaxy 
jPuerarilll993l) . this corresponds to a ~ 6° azimuthal offset at a 



Soft X-ray emission from the inner disk ofM33 1 1 
M33 



O 



o - 



is 




-t — I — I — I — I — I — I — I — I — I — I — I — I — I — I — h 




H 1 1 1 1 1 1 1 1 1 1 1 h 



"L+j-i-h— K+" 
k ' ' ' ' k 



-l""Li_rH — H 



_i i i i i i i i i i i i i i i i_ 



Radial distance (kpc) 

Figure 4. A comparison of the radial profile of X-ray emission with FUV emission, SFR and mass distributions in the central disk of M33. The x-axis 
refers to the major-axis radius scaled to kpc, assuming the distance to M33 is 795 kpc. The following information is provided: Top panel: The soft X-ray 
surface brightness versus radius (upper curve). The radial profile of the FUV emission from GALEX images (middle curve). The radial profile of the K-band 
emission from 2MASS images (lower curve). Second panel: The ratio of the soft X-ray luminosity in erg s — 1 pc — 2 (0.3-1 keV) to the local SFR in units of 
Mq yr — 1 pc -2 . Third panel: The ratio of the soft X-ray luminosity in erg s — 1 pc~ 2 (0.3-1 keV) to the K-band derived mass in units of Mq pc -2 . Bottom 
panel: Variation in X-ray spectral hardness, (H-S)/(H+S), versus radius, where H refers to the 0.8—1.2 keV band and S to the 0.3-0.8 keV band. 



galactocentric radius of 2 kpc, between the X-ray emission and the 
FUV emission (with the former leading the latter since at 2 kpc 
we are within the corotational radius). Such an offset would not 
be observable in M33 in XMM-Newton data and, in fact, such ef- 
fects have yet to be identified even in the hig her spatial resolut ion 
Chandra observations of other spiral galaxies dTvler et alj2004T) . 

Even after allowing for the contribution of the old stellar popu- 
lation, we are left with a substantial fraction of the soft X-ray emis- 
sion originating in a smoothly distributed component. Since M33 
has an intermediate inclination, it is not possible to disentangle disk 
emissi on from a lower-halo component. However. IStrickland et al.l 
(2000), in a study of 9 nearby edge-on spiral starburst and normal 
galaxies, have managed to distinguish between unresolved disk and 



lower-halo components and find that the former dominates across 
their sample. In all but the systems with highest SFR, the luminos- 
ity of the halo emission is at least a factor of 3-4 lower than that 
observed in the disk. In the case of M33 it therefore seems reason- 
able to assume that the bulk of the smoothly-varying component is 
confined to the disk but located both in the spiral arms and inter-arm 
regions. A cooling timescale of ~ 10 9 7? 1//2 yr for this component 
is consistent with its smooth distribution across the inner disk of 
M33. 
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7 CONCLUSIONS 

We have used archival XMM-Newton observations to examine the 
residual X-ray emission observed from M33 after the exclusion 
of the bright point source population to a limit of Lx > 2 x 
10 35 erg s _1 . Using the same methodology as in OW09, we have 
investigated the spectral and spatial properties of the X-ray emis- 
sion within an inner disk region extending up to 3.5 kpc from the 
nucleus of the galaxy. 

The observed X-ray spectrum can be modelled as thermal 
emission with cool and hot components of ~ 0.2 keV and m 0.6 
keV respectively, with the cooler emission providing the dominant 
contribution to the luminosity. There is some evidence for a subso- 
lar metallicity consistent with other indicators of low-metallicity in 
M33. 

The strong correlation established between X-ray and FUV 
morphologies confirms the close linkage between X-ray emis- 
sion and recent star formation. Detailed comparison of soft X- 
ray and FUV radial profiles in the inner disk of M33 reveals the 
ratio of X-ray emission to SFR to be 1 — 1.5 x 10 39 erg s _1 
( Mq yr~ )~ (in the 0.3-2 keV band). This match es the pre- 
dictions of iMas-Hesse. Oti-Floranes~& Cervind (2008) for an ex- 
tended burst of star formation occurring 10 Myr ago, with an ef- 
ficiency of mechanical energy conversion to X-rays of ~ 1%. 
The soft X-ray emission to mass ratio found for M33 is 4 x 
10 2S erg s _1 M0 -1 , a factor of 5 higher than the corresponding 
value for dwarf elliptical galaxies and spirals with low SFR. This 
implies that up to ~ 20% of the observed soft X-ray emission orig- 
inates in the old stellar source population, in source types such as 
cataclysmic variables and active binaries. 

With the contribution of the old stellar population subtracted, 
the soft X-ray emission is found to be equally split between two 
spatial components, one which closely traces the spiral arms of the 
galaxy and the other more smoothly distributed across the inner 
disk region. The constraints on the cooling timescale implied by 
the presence of spiral features in the soft X-ray images suggest the 
presence of a highly clumped component, encompassing sites of 
on-going starformation, HII regions and hot gas bubbles. The nature 
of the smoothly distributed component is much less certain. Plausi- 
bly it may represent the integrated emission from a whole range of 
processes and source types including supernovae occurring in the 
interarm regions, individual sources and source complexes with lu- 
minosity not far below the applied luminosity threshold and an ac- 
cumulated distribution of hot gas which has managed to leak away 
from the site of its original production. 
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